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Abstract — We theoretically investigate the 
recent photoemission-type experiment on 40 K Fermi gases 
done by JILA group. Including pairing fluctuations within 
a strong-coupling T-matrix theory, as well as effects of a 
harmonic trap within the local density approximation, we 
calculate photoemission spectra in the BCS (Bardeen-Cooper- 
Schrieffer)-BEC (Bose-Einstein condensation) crossover re- 
gion. We show that the energy resolution of the current 
photoemission experiment is enough to detect the pseudogap 
phenomenon. We also show how the pseudogap in single- 
particle excitations continuously changes into the superfluid 
gap, as one decreases the temperature below the superfluid 
phase transition temperature. Our results would be useful for 
the study of single-particle properties of ultracold Fermi gases 
in the BCS-BEC crossover. 

Index Terms — Atomic Fermi Gas, Pseudogap, Superfluidity 



I. Introduction 

THE recent photoemission-type experiment done by 
JILA group[l], [2] is a very powerful method to exam- 
ine single-particle properties of ultracold Fermi gases. This 
Fermi system has the unique property that the strength of a 
pairing interaction can be tuned by adjusting the threshold 
energy of a Feshbach resonance[3], [4]. Using this, one can 
study from the weak-coupling Bardeen-Cooper-Schrieffer 
(BCS)-type superfluid to the Bose-Einstein condensation 
(BEC) of tightly bound molecules in a unified manner[3], 
[4], [5], [6], [7], [8], [9], [10], [11]. The continuous 
change from the weak-coupling BCS regime to the strong- 
coupling BEC regime is frequently referred to as the BCS- 
BEC crossover phenomenon, which is one of the most 
exciting topics in cold atom physics [3], [4]. In particular, 
the so-called pseudogap phenomenon has recently attracted 
much attention[12], [13], [14], [15], [16], [17], [18], [19], 
[20], because this strong-coupling phenomenon is also 
considered as a key to clarify the mechanism of high-T c 
cuprates[21]. Since the pseudogap appears in single-particle 
excitation spectra, the photoemission-type experiment de- 
veloped by JILA group[l], [2] is very suitable for the study 
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of this problem. Indeed, the observed excitation spectra in 
the crossover regime clearly deviate from those in the case 
of a free Fermi gas[l], [2]. 

In high-T c cuprates, the mechanism of pseudogap has not 
been completely clarified yet, because of the complexity of 
the system (although the importance of pairing fluctuations 
has been pointed out[21]). On the other hand, the origin of 
the pseudogap is well known in cold Fermi gases, namely, 
the formation of preformed Cooper pairs by strong pairing 
interaction. Thus, in the latter system, one can conveniently 
examine pseudogap physics originating from superfluid 
fluctuations, without any ambiguity. 

In considering the photoemission-type experiment on 
cold Fermi gases, one should note the following two ex- 
perimental situations. First, the system is always trapped in 
a harmonic potential [16], [20]. As a result, pairing fluctua- 
tions are spatially inhomogeneous, which naturally leads to 
inhomogeneous pseudogap effects. For example, one may 
expect the situation that, while the pseudogap is remarkable 
in the trap center, such an effect is weak around the edge 
of the trap. Since the current photoemission experiment 
does not have spatial resolution[l], it only gives spatially 
averaged data. Thus, it is an interesting problem how the 
pseudogap effect can be seen in such spatially averaged 
spectra. Second, the observed photoemission spectrum is 
always affected by finite energy resolution[16]. Because of 
this, the pseudogap structure in the spectrum is expected to 
be smeared to some extent. Thus, it is an important problem 
whether or not the current experiment with finite energy 
resolution (D ~ 0.2£p, where ep is the Fermi energy) is 
enough to detect the interesting pseudogap phenomenon. 

Between the above two important issues, we have exam- 
ined the first one in a previous paper[17]. Including pairing 
fluctuations within a strong-coupling T-matrix theory, as 
well as effects of a harmonic trap within the local density 
approximation (LDA), we showed how the pseudogap 
appears in the spatially averaged photoemission spectra at 
the superfluid transition temperature T c . In this paper, we 
extend our previous work[17] to include the experimental 
energy resolution. We show that the recent photoemission 
measurements on 40 K Fermi gases[l], [2] really detect the 
pseudogap. We further extend this work to the superfluid 
state below T c . We clarify how the pseudogap continuously 
changes into the superfluid gap below T c . Far below T c , the 
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single-particle spectral weight is shown to exhibit a double 
peak structure, originating from quantum fluctuations and 
superfluid gap. 

The outline of this paper is as follows. In Sec. II, we 
present our formulation. We explain how to include strong 
pairing fluctuations, effects of a harmonic trap, as well as 
finite energy resolution. In Sec. Ill, we show our numerical 
results on photoemission spectra. Throughout this paper, 
we take h = k& = 1. 

II. Formalism 

We consider a two-component Femi gas, described by 
pseudo spin a =t, I- In real cold Fermi gases, these pseu- 
dospins physically represent two atomic hyperfine states 
contributing to the formation of Cooper pairs. So far, all 
the current experiments on cold Fermi gases are using a 
broad Feshbach resonance to tune the strength of a pairing 
interaction[l], [2]. In this case, it is known that the pairing 
interaction may be simply treated by the ordinary BCS 
model[3], as far as we consider the interesting BCS-BEC 
crossover physics. The model Hamiltonian is given by [17], 
[18], [19], [23] 

H = E*^ 7 * - Ar ^P -U^p + (q)p-(q). (1) 

P Q 

Here, ^ P = (c^,c_ p ^) is the two-component Nambu 
field, where c pa is the annihilation operator of a Fermi atom 
with pseudospin a. Tj (j = 1, 2, 3) the Pauli matrices acting 
on the particle-hole space. £ p = e p — p is the kinetic energy 
e p of a Fermi atom, measured from the chemical potential 
p. The pairing interaction is described by —U, which is 
assumed to be a tunable parameter. In cold atom physics, 
the strength of an interaction is usually measured in terms 
of the observable scattering length a s . In the present BCS 
model, it is related to — U as [22], 
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where m is an atomic mass, and uj c is a high-energy 
cutoff. In this scale, the weak-coupling BCS regime and 
the strong-coupling BEC regime are, respectively, given 

by (/cf^s) -1 ^ — 1 an d (&Fa s ) -1 ~ 1 (where fcp is 
the Fermi momentum). The intermediate coupling region, 
— 1 < (&Fa s ) -1 ~ 1, is the crossover region. 

The superfluid phase is characterized by the superfluid 
order parameter, A = ^ p ( c -?4 c pt) (which is taken to be 
real, and is proportional to the T\ -component in Eq. (1). 
In the last term of Eq.(l), p±(q) = [pi(q) ± ip2(q)]/2 in 
Eq. (1) involves the generalized density operators Pj(q) = 
^p+g/2 r ^P-<?/2 0' = 1,2), describing amplitude 
(j = 1) and phase (j = 2) fluctuations of the order 
parameter. 

We treat the interaction term — U in Eq. (1) within the T- 
matrix approximation [12], [13], [14], [17], [18], [19]. For 
this purpose, we introduce the 2 x 2-matrix single-particle 



thermal Green's function, 



Gp(iuj n ) 1 - Yip(iuj n ,r)' 
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Here, Gp(iuj n )~ l = iuj n — + Ati is the mean- 
field Green's function, where uj n is the Fermi Matsubara 
frequency. The 2 x 2-matrix self-energy Yip(iuj n ), involving 
effects of pairing fluctuations within the T-matrix approx- 
imation, has the form[19] 

Yip{iuj n ) 

= E K S \^n)r-sGl +q {iu n ^iv n )r^, 



q,v n s,s f 



(4) 



where r± = t\ ± iT2, and v n is the Bose Matsubara 
frequency. The scattering matrix in the Cooper channel 
T ss ' (iv n ) is given by[19] 
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Here, 



(6) 

is the lowest order of the pair-correlation function in terms 
of the interaction —U. 

Experimentally, since a Fermi gas is always trapped in 
a harmonic potential V(r) = rac^ r r 2 /2[24], we include 
this inhomogeneous effect within the local density approx- 
imation (LDA), which is simply achieved by replacing the 
chemical potential p by p(r) = p — VtrapM- The single- 
particle Green's function and superfluid order parameter 
then have spatial dependences as G p (iuj n ,r) and A(r), 
respectively. 

The LDA photoemission spectrum I(p, fi) in a trapped 
Fermi gas is given by [18] 



J(p,fi) 
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'drp 2 Ap(t P (r)-n,r)f(Z p (r)-n), 

(7) 

where £ p (r) = e p — p(r), f(oj) is the Fermi distribution 
function, and a is a constant factor. (The detailed expres- 
sion of a is not necessary in the following discussions.) 
The single-particle spectral weight A p (^ p {r) — f£, r) at r is 
obtained from the analytic continuation of the LDA Green's 
function as 



A p (u,r) = — — lmG p (iuj ri 

7T 



=u + i6,r)\u. (8) 



As mentioned previously, the observed photoemission 
spectra in Ref. [1] are affected by experimental energy 
resolution. Incorporating this situation into Eq. (8), we have 

J(p,fi) = / dzI(p,Sl)e-^-. (9) 

V 71 "^ J-00 
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For the value of energy resolution D, we employ the 
experimental value D ~ 0.2£p[l]. Equation (9) is di- 
rectly related to the occupied single-particle spectral weight 
5(p,cj)as[18] 



S(p,u>) = -I(p,Sl- 
a 



(10) 



In the simplest non-interacting Fermi gas at T = 0, 
S(p,uj) has a J-functional peak line along the free-particle 
dispersion uj = £ p when uo < 0. The spectral weight in the 
positive energy region vanishes due to the vanishing Fermi 
distribution function f(u > 0) = at T = 0. 

To calculate Eq. (9), we first determine the superfluid 
order parameter A(r) and chemical potential /i, which is 
achieved by solving the gap equation 



1 = 
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/ tanhjgj 



= 0, (11) 



(where E p (r) = v / '^ P (r) 2 + A(r) 2 is the Bogoliubov 
single-particle excitation spectrum), together with the LDA 
number equation for Fermi atoms (within the T-matrix 
approximation), 
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The superfluid phase transition temperature T c is obtained 
under the condition that the gap equation (11) is satisfied 
when A(r = 0) = 0. Using the self-consistent solutions, 
we calculate the spectral weight in Eq. (10). 

III. Photoemission spectra with finite energy 

RESOLUTION IN THE BCS-BEC CROSSOVER 

Figure 1 shows the occupied single-particle spectral 
weight S(p,uj) at T c . Comparing this figure with the case 
with no energy resolution (D = 0) shown in Fig. 2, we 
find that the finite energy resolution (D = 0.2ev) broadens 
the spectra. However, even in this case, one can still see 
the pseudogap effect in Fig.l. That is, starting from the 
weak-coupling BCS regime shown in Fig. 1(a), we see that 
the spectral peak gradually deviates from the free particle 
dispersion, as the pairing interaction becomes strong. (See 
Fig. 1(b).) In the strong-coupling BEC regime (Fig. 1(c)), 
the single-particle excitation spectrum splits into an upper 
sharp particle branch (uo > 0) and lower broad hole branch 
(oj < 0). Since the superfluid order parameter is absent 
at T c , the pseudogap structure in Fig. 1(c) purely comes 
from strong pairing fluctuations. As shown in Fig. 1(b), the 
calculated lower peak line (white dashed line in the figure) 
agrees well with the recent experiment on a 40 K Fermi 
gas done by JILA group. As discussed in our previous 
papers [17], [18], [19], this lower peak is characteristic 
of the pseudogap effect, originating from a particle-hole 
coupling induced by pairing fluctuations. (Note that, when 
the pairing interaction is absent, the photoemission spec- 
trum only has the upper peak (black solid line in Fig. 
1(b)) line along the dispersion of a free atoms. These 




Fig. 1. Calculated intensity of occupied spectral weight S(p, ou) in a 
trapped Fermi gas at T c . The bright color region shows high intensity. 
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panel (b), the black solid line shows the upper peak position of the 
spectrum, and white dashed line shows the lower peak positions of the 
spectrum. (Although the peak intensity of the lower line cannot be seen in 
panel (b), this is simply because its magnitude is much smaller than the 
dominant peak intensity around p/k F ~ 0.5.) The white dots represent 
the experimental data for the lower peak positions measured in Ref.[l]. 



results naturally lead to the conclusion that the energy 
resolution D ~ 0.2sf at the current stage of photoemission 
experiment is enough to detect the pseudogap phenomena 
in the BCS-BEC crossover regime of cold Fermi gases. 

We now proceed to the superfluid phase below T c . Far 
below T c , since thermal pairing fluctuations are almost 
absent, the occupied spectral weight S(p,uj) in the positive 
energy region almost vanishes, as shown in Fig. 3. In the 
ordinary mean- field BCS theory for a uniform Fermi super- 
fluid, the single-particle Bogoliubov excitation spectrum in 
the hole branch is given by 



-E n 



e P 



A 2 . 



(13) 



Although thermal fluctuations are almost absent far below 
T c , Fig. 3 shows that the peak positions of the occupied 
spectral weight still deviate from the expected mean-field 
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result in Eq. (13). In particular, one sees two branches in 
panels (b) and (c). 

To see the origin of the appearance of two branches 
in Figs. 3(b) and (c), we show in Fig. 4 the temperature 
dependence of the occupied spectral weight in the crossover 
region (l/k^a s = 0.15). As one decreases the temperature 
below T c , panel (b) indicates that the lower peak line 
(dashed line) is pushed down by the development of 
superfluid order parameter. At lower temperatures shown 
in panel (c), the other peak line (green dotted line) ap- 
pears, which gradually reduces to the dispersion of hole 
Bogoliubov excitations in Eq. (13), as shown in panels (d) 
and (e). That is, between the two branches in Figs. 3(b) and 
(c), the lower one is related to the pseudogap at T c , and 
the upper one originates from the ordinary BCS excitation 
gap. As mentioned previously, since thermal fluctuations 
are almost absent far below T c , the lower peak at T <C T c 
is considered to be also related to quantum fluctuations. 

In a previous paper[18], we showed that, in the BEC 
regime, the lower peak energies in S(p,u) calculated at 
T c are slightly larger than the experimental data by JILA 
group [1]. In this regard, we note that it is difficult to 



accurately determine the temperature in cold Fermi gases. 
Thus, since the lower peak line is pushed down in the 
superfluid phase (See Fig. 4.), a possible idea to resolve this 
discrepancy is that the experiment in the BEC regime was 
actually done in the superfluid phase below T c . Indeed, 
when we take T — 0.6T C < T c , our theoretical result 
in the BEC regime (l/(k^a s ) = +1) well explains the 
observed peak energies, as shown in Fig. 5. Although we 
need further studies about this problem, our result shows 
that the existence of a finite superfluid order parameter 
is a possible idea to explain the photoemission spectrum 
observed in the BEC regime. 

To summarize the behavior of S(p,u) in the BCS- 
BEC crossover below T c , we introduce three characteristic 
temperatures shown in Fig. 6. Just below T c , although the 
system is in the superfluid phase, one cannot see any su- 
perfluid effect in S(p,u). Below a certain temperature Ti, 
a peak structure corresponding to the BCS excitation gap 
appears in the spectrum. One may also introduce another 
characteristic temperature T2, below which this peak line 
is well described by the Bogoliubov excitation spectrum 
in Eq. (13). Since thermal excitations are suppressed far 




Fig. 5. Calculated occupied single-particle spectral weight S(p, lj) in 
the BEC regime at (k^ds)' 1 = 1. We take at T = 0.6T C . For clarity, 
we also draw three lines at the peak positions of the spectral weight. The 
white dots are experimental data at (fepOs) -1 = 1[1]. 



Fig. 4. Calculated temperature dependence of occupied single-particle 
spectral weight S(p, lj). We set l/(k-pa s ) = 0.15. The lines shown 
in this figure show peak positions of the spectral weight. The calculated 
chemical potential /i are (a)0.33e F , (b)0.43e F , (c)0.47e F , (d)0.56e F , and 
(e)0.60e F . 
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Fig. 6. Characteristic temperatures obtained from the occupied single- 
particle spectral weight S(p, lj). Below T\, the peak line associated with 
the BCS excitation gap appears in S(p, lj). The left edge of this peak line 
reaches p = at T2. Below T2, this line looks similar to the dispersion 
of hole Bogoliubov excitations in Eq. (13). Below T3, the spectral weight 
no longer has a finite weight in the positive energy region within our 
numerical accuracy. In this figure, we also plot T = 2|/x(T c )| when 
/i < 0. Since 2|/x(T c )| gives the binding energy of a two-body bound 
state, it physically means the characteristic temperature, above which two- 
bound states dissociate into atoms thermally. Thus, the right side of this 
line may be viewed as the region of a gas of tightly bound molecules, 
rather than strongly correlated Fermi atoms. 



below T c , S(p J uj) becomes absent in the positive energy 
region u > 0. The characteristic temperature T3 is defined 
as the temperature when this situation is realized. 

In Fig. 6, we also plot the binding energy of a two- 
body bound molecule, give by E g = 2|/i(T c )|[17], [18], 
[19]. When the temperature is lower than E g , thermal 
dissociations of molecules are suppressed, so that the 
system is close to a molecular Bose gas formed by a two- 
body effect, rather than a strongly correlated Fermi gas. 

IV. Summary 

To summarize, we have theoretically investigated the 
recently realized photoemission-type experiment on ultra- 
cold Fermi gases. Our theory takes into account pairing 
fluctuations within the T-matrix, as well as effects of a 
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harmonic trap within the local density approximation. To 
include the experimental situation, we have also taken into 
account the energy resolution D = 0.2£p in calculating the 
photoemission spectrum. 

At T c , we showed that, although the energy resolution 
broadens the photoemission spectra, we can still see the 
pseudogap effect in the BCS-BEC crossover region. Below 
T c , we also showed how the pseudogapped spectral weight 
continuously changes into the spectral weight with a finite 
BCS superfluid gap. To summarize the temperature depen- 
dence of spectral weight below T c , we have introduced 
three characteristic temperatures, Ti, T2 and T3. Since the 
photoemission-type experiment is one of the most powerful 
methods to observe single-particle properties of cold Fermi 
gases, our results would be useful for the understanding of 
strong-coupling effects in the BCS-BEC crossover regime 
of this system. 
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